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ABSTRACT 

The  sound  design  of  environmental  control  systems  lor  underground 
shelters  is  based  on  knowledge  of  the  thermal  transient  response  under 
certain  climatic  conditions.  To  gain  insight  into  the  transient 
response  of  underground  shelters,  a  series  of  tests  using scale 
model  shelter  was  conducted.  Model-prototype  considerations  resulted 
in  several  model  temperature  distortions  which  could  be  accounted  for  by 
analytical  techniques.  The  corrected  model  results  agreed  well  with 
results  from  an  analog  computer  study  which  considered  the  same  prototype 
shelter.  The  model  results  also  Indicated  that  shelter  shape  does  not 
significantly  affect  its  transient  response. 

Another  phase  of  this  study  was  the  development  of  a  non-computer 
design  procedure  for  determining  the  environmental  control  system 
capacity  required  for  a  given  set  of  climatic  and  soil  conditions.  The 
design  procedure  was  applied  to  a  number  of  shelter  locations  and  cli¬ 
matic  conditions  to  test  its  performance.  As  expected,  the  solutions 
indicated  that  ventilation  rate  and  air-conditioning  capacity  depend 
heavily  on  climate,  initial  soil  temperature,  and  shelter  area  per 
person. 

This  report  is  Intended  as  a  record  of  the  analytical  and  experimental 
methods  evolved  under  the  task  to  date.  All  analytical  and  experimental 
techniques  are  described  in  detail,  and  numerical  examples  for  a  typical 
model  run  and  design  procedure  solution  are  given.  Further  refinement 
and  application  of  rhe  model  study  techniques  and  non-computer  design 
procedure  are  presently  underway,  and  it  is  hoped  that  this  work  will 
lead  to  recommended  design  data. 


Qualified  requesters  may  obtain  copies  of  this  document  from  DDC. 
Release  to  the  Clearinghouse  is  Authorized. 


NOMENCLATURE 


Primary  symbols  (presented  in  order  of  appearance  in  text). 

x,y,z  *  distance  into  the  thermal  medium  (soil);  ft. 
q  »  heat  rate;  Btu/hr,  Btu/ hr -person, 
t  ■  temperatures;  °F 

k  =  thermal  conductivity;  Btu/hr-ft-°F. 
p  ■  density;  Ib/ft^. 
c  =  specific  heat;  Btu/lb-°F. 

Q  =  time;  hours,  days. 

A  »  area;  ft  ,  ft  /person. 

V  ■  volume;  ft^. 

L  *  any  length  corresponding  on  the  model  and  prototype;  ft. 
OL  *  thermal  diffusivity;  ft  /hr. 

N  »  primary  non-dimensional  parameter. 

h  “  average  convective  heat  transfer  coefficient  along  the 
boundary  layer;  Btu/hr-ft^-°F. 

Bi  =  Biot  number;  hL/k^,  (non-dimensional). 

O 

Q  ■  ventilation  rate;  ft  /hr-person. 

R  *  air-conditioning  capacity;  tons/person,  Btu/hr-person. 

F  *  sensible  heat  factor;  non-dimensional, 
w  *  absolute  humidity;  lb  moisture/lb  dry  air. 

H  ■  enthalpy  of  condensation;  Btu/lb. 
n  ■  number  of  shelter  occupants. 


Hi 


Superscripts 


*  denotes  a  non-dimensional  ratio;  ()/(). 

m  p 

Subscripts 

x.y.z  *  quantities  with  respect  to  distance  in  the  thermal  medium, 
s  *  a  stored  quantity. 

w  »  quantities  with  respect  to  the  shelter  wall  (soil), 

m  =  quantities  with  respect  to  the  model, 

p  =  quantities  with  respect  to  the  prototype, 

a  *  quantities  with  respect  to  air  inside  the  shelter. 

A  >  quantities  with  respect  to  air  outside  the  shelter, 

o  ■  an  initial  condition;  6  ■  0. 

0  ■  a  condition  at  time  6*^0. 

S  >  quantities  with  respect  to  sensible  heat. 

L  ■  quantities  with  respect  to  latent  heat. 

I  ■  a  total  quantity. 

g  >  quantities  with  respect  to  heat  generated, 
r  ■  quantities  with  respect  to  air  conditioning. 

V  -  quantities  with  respect  to  ventilation. 


NOTE:  All  other  symbols  are  defined  as  they  appear  in  the  text. 
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INTRODUCTION 


Underground  blast  and  fallout  shelters  are  usually  designed  to 
pr(.'vitlc  prot^rMon  and  habitable  living  conditions  at  the  lowest 
reasonable  cost  per  person.  The  requirement  for  economical  anderQioa.^H 
shelter  facilities  has  led  to  recognizing  the  soil  surrounding  the  shel¬ 
ter  as  a  heat  sink  which  can  allow  a  significant  reduction  In  the 
ventilation  and  air-conditloning  requirements.  The  determination  of  the 
required  environmental  control  system  capacity,  Including  heat  absorption 
by  the  soil,  Is  a  difficult  problem,  but  digital  and  analog  computer 
solutions  have  been  devised  to  solve  lt.^>^  These  solutions  give  the 
thermal  transient  response  of  shelters  under  certain  operating  and 
climatic  conditions.  Environmental  control  studies  have  also  been  con¬ 
ducted  with  actual  underground  shelters. Unfortunately  tests  with 
actual  shelters  have  been  limited.  All  of  these  studies  (Including 
similar  work  not  referenced)  has  given  considerable  Insight  Into  the 
environmental  control  of  underground  shelters;  however  it  was  felt  that 
additional  Insight  and  evidence  was  needed  in  view  of  the  limited 
experience  with  actual  shelters.  It  was  noted  that  this  Information, 
coupled  with  the  results  of  prior  work,  might  provide  the  basis  for  a 
non-computer  procedure  for  the  design  of  underground  shelter  environmental 
control  systems  Since  the  majority  of  prior  work  was  basically  analytical, 
it  was  decided  that  additional  analytical  work  would  be  duplication  of 
effort.  On  the  other  hand,  tests  with  actual  shelters  are  expensive  and 
unwieldy  experimentally  and  loglstlcally .  Thus  It  was  decided  to  approach 
the  shelter  problem  by  means  of  model  studies.  The  model  would  provide 
a  degree  of  realism  which  cannot  be  easily  achieved  with  analySi.s.  For 
example,  the  shape  of  the  shelter  could  be  duplicated,  and  heat  flow 
through  the  soil  would  actually  occur  In  three  dimensions  (most  of  the 
prior  analytical  work  assumed  an  idealized  shelter  shape  and  heat  flow  In 
one  dimension). 

To  simplify  construction  of  the  model  and  reduce  experimental 
complexities.  It  was  decided  to  model  a  steel  arch  structure.  This  type 
of  shelter  has  been  considered  in  previous  work,  and  It  does  not 

involve  consideration  of  large  concrete  masses  relative  to  soil  Msses. 

The  model  shelter  was  burled  In  a  suitable  soil  medium  and  heated 
appropriately  to  simulate  ventilation,  alr-conditloning,  and  the  physio¬ 
logical  thermal  response  to  dry-bulb  temperature  of  people  engaged  In 
moderate  activity.  Temperatures  were  measured  by  thermocouples  placed 
Inside  the  shelter,  on  the  shelter  walls,  and  In  the  soil  surrounding 
the  shelter.  The  model  exhibited  temperature  distortions  because  of  the 
method  of  heat  Input,  and  because  not  all  of  the  non-dlownslonal  parameters 


could  be  adjusted  to  satisfy  the  laws  of  similitude.  However,  all  of 
these  distortions  were  accounted  for  analytically.  The  results  of  the 
model  experiments,  when  distortions  were  corrected,  showed  excellent 
agreement  with  the  analog  computer  solutions. ^  A  quantitative  comparison 
of  the  model  behavior  with  results  from  other  work  could  not  be  performed 
because  of  insufficient  data  reported;  nevertheless  the  results  were 
similar  in  general  terms.  In  general,  the  model  showed  th^lt:  (1)  after 
occupation  the  air  temperature  rises  very  sharply  during  the  first  day 
after  which  there  is  a  marked  decrease  in  temperature  rise  and  the  response 
is  quite  flat  --  approximating  a  step  function  --  however  a  steady-state 
condition  is  not  usually  reached  during  the  expected  period  of  occupancy 
(10  to  14  days);  (2)  the  shape  of  the  shelter  does  not  greatly  affect  its 
transient  response,  thus  prior  assumptions  of  one-dimensional  heat  flow 
appear  to  be  reasonable.  An  unanticipated  result  of  the  model  study  is 
that  the  distortion  analysis  can  be  applied  to  full-size  shelters,  which 
raises  the  possibility  that  a  program  of  environmental  control  studies 
covering  a  range  of  locations  and  climatic  conditions  might  be  undertaken 
on  a  single  prototype  shelter  installation. 

After  the  model  studies  indicated  that  the  geometrical  shape  of  the 
shelter  does  not  significantly  affect  its  transient  response,  and  that 
the  response  might  be  approximated  by  a  step  function,  it  was  found  that 
a  non-computer  procedure  for  the  design  of  underground  shelter  environ¬ 
mental  control  systems  could  be  developed.  The  procedure  involves  a 
certain  degree  of  trial  and  error  and  a  working  knowledge  of  paychro- 
metric  charts;  otherwise  it  is  conceptually  straightforward,  flexible, 
and  can  be  used  without  difficulty  by  trained  personnel.  During  the  early 
portion  of  the  transient  response,  the  design  procedure  gives  only  fairly 
accurate  results  because  the  assumption  of  a  step  function  response  is 
somewhat  in  error.  However,  the  time  of  interest  is  really  at  the  end  of 
the  expected  occupancy  period,  and  at  this  point  the  procedure  usually 
gives  results  within  l^F  of  the  analog  computer  and  model  stuvdy  results. 

This  report  considers  the  experimental  and  analytical  methods  used 
to  develop  the  model  study  and  non-computer  design  procedure,  and  it  is 
divided  into  three  parts.  The  first  part  deals  with  the  model  study,  the 
second  part  is  devoted  to  the  non-computer  design  procedure,  and  the  third 
part  gives  general  recommendations.  The  design  procedure  is  applied  to  a 
number  of  shelter  locations  and  climatic  conditions;  however  the  results 
are  to  be  used  only  for  purposes  of  comparison.  Further  refinement  and 
application  of  the  design  procedure  must  be  undertaken  before  it  can  be 
used  with  confidence.  This  work  is  receiving  continued  attention,  and  It 
will  be  reported  at  a  latnir  date. 
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PART  I.  THE  MODEL  STUDY 


Development  of  the  Non-Dimensional  Parameters 

The  non-dimensional  parameters  used  for  the  shelter  model  study 
can  be  derived  conveniently  by  considering  the  heat  transfer  phenomena 
involved.  The  nomenclature  used  is  listed  at  the  end  of  the  text,  and 
the  assumptions  taken  are  given  in  Appendix  A. 


To  begin  the  analysis,  consider  a  differential  element  in  cartesian 
coordinates  with  heat  flow  in  all  three  dimensions  (see  Figure  1).  The 
energy  balance  between  _ lactsd  end  heat  stored  is: 

(d,^  +  d,y  +  dq^)  -  dq^  +  (dq_^^_^  +  d^^j^  +  dq^_^^^)  (1 


At  the  surfaces  of  the  differential  volume,  the  Fourier  conduction 
equation  is  expressed  in  terms  of  partial  differential  equations: 

■'"x  ■  It 


dq 


x+dx 


■k  (dy)  (da) 
w 


(2) 

(3) 


(Similar  relations  can  be  developed  for  dq^,  ^**rfda^ 

The  heat  stored  within  the  volume  element,  dq  ,  is  expressed  as 

8 


dq 


8 


(^) 


To  expedite  the  analysis,  let  dA  ■  (dx)  (dy)  ■  (dy)  (da)  -  (dx)  (dc)  and 
dV  -  (dx)  (dy)  (da). 


Now,  by  substitution.  Equation  1  can  be  expressed  as 


-k„«iA) 


^  dy^  da^ 


(5) 


3 


t 


The  first  order  terms  containing*^,  and  can  be  eliminated 
immediately  so  that 

^^t  ,  ^^t  ,  ^^t  ,  . 

— r  dx  + - r  dy  + - r  dz 

dz 


k  (dA) 
w 


(6) 


Equation  6  is  dimensionally  consistent;  therefore  the  quantity,  L,  can  be 
substituted  for  all  quantities  involving  length,  and  Equation  6  can  be 
expressed  dimensionally  as 


L 


P  c  L' 

w  w 


9 


(7) 


which  can  be  simplified  to 


P  c  L 
w  w 

e 


(7a) 


It  should  be  noted  that  the  quantities  k^,  p^,  and  c^  have  not  been 

expressed  in  dimensional  form  because  these  quantities  are  recut^nizcd 
soil  properties.  Equation  7a  can  be  arranged  to  yield  a  nondimens ional 
parameter. 


a  9 


where  a  •  (■~) 

w  pc  w 


A  basic  principle  of  similitude  is  that 


^  ^  O'  ^ 

) 

L  model  L  prototype 


(8) 


(9) 


The  quantity,  L,  is  a  convenient  length  for  proportionately  scaled  nx>dels 
It  does  not  matter  which  length  is  chosen,  but  the  length  chosen  on  the 
prototype  must  correspond  to  the  length  taken  on  the  model. 


The  parameter,  N.,  has  been  derived  from  the  interaction  between  heat 
conducted  and  halt  stored.  It  is  this  Interaction  which  affects  the 
transient  temperature  response  in  the  thermal  medium.  N  can  be  used  to 
determine  prototype  time  values  for  corresponding  model  ^time  values. 
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Let  L*  -  L  /L 
m  p 


and  a  *  “  a  fit 
w  w  /  w 
m  p 


Equation  9  can  be  solved  for  0  . 

P 


6 

P 


a  * 


(10) 


A  non-dimensional  parameter  for  temperatures  in  the  soil  medium  at 
any  point  can  be  derived  from  the  Fourier  conduction  equation. 


dq 


-k  (dA) 
w 


dt  dt  dt 

_ w  ,  w  w 

sr 


(11) 


e  quantities  q,  k  , 
for  A,  and  L  for  x, 


and  t  will  be  retained,  and  upon  substitution  of 
y,  and  z,  Equation  11  can  be  expressed  dimensionally 


as 


q 


d  2 
•  k  L 
w 


C 

ja 

L 


The  resulting  non-dimensional  parameter  is 


(12) 


N 


7 


k  Lt 


w _ y 


q 


(13) 


The  temperature,  t^,  varies  as  a  function  of  time,  6,  from  some  datum 

temperature,  t^  ,  at  9  ■  0.  Therefore,  can  be  expressed  more 
o 

conveniently  as 


k  L(t  -t  ) 

WWW 

"2  ■ - T— ^ 


Applying  the  basic  principles  of  similitude. 


k  L(t  -t  ) 
w  w  w 

_ a. 


model 


K  L(t  -t  ) 
w  w  w 
_ 2_ 


prototype 


(I3a) 


(U) 


5 


If  we  let  q*  =  q  /q  and  k  *  =  k  /k  ,  the  prototype  shelter  wail 

^  m  p 

temperature  rise,  (t  -t  ),  can  be  expressed  as 

o 


k  *L* 

(t  -t  )  =  —  (t  -t  )m 

ww  p  q*  ww 


(15) 


Heat  tl')w  from  the  shelter  air  to  the  shelter  wall  is  expressed  by 
the  familia*.  convection  “qu  ^on  written  below  in  differential  form. 


dq  =  h(dA)  (t^-t^) 


(16) 


where:  h  is  the  average  convective  heat  transfer  coefficient  along 
the  boundary  layer. 

t  is  the  air  temperature  beyond  the  boundary  layer. 

r 

This  equation  immediately  suggests  a  non-dimensional  parameter  when  L 
is  substituted  for  A. 


N-  =  —  (t  -t  ) 
3  q  a  w 


(17) 


Le! ting  h*  =  h  /h  and  imposing  the  principle  of  similitude,  we  have  an 
m  p 

expression  for  the  prototype  air-wall  temperature  difference. 

,2 


(t  -t )  -t ) 

a  w  p  q*  a  w  m 


n 


Model-Prototype  Temperature  Distortions 

Equations  15  and  18  are  basic  temperature  relations  between  the 
prototype  and  model.  I'Jhen  considered  together,  these  equations  yield  an 
important  result.  Recall  Equation  15  and  compare  it  to  Equation  18.  Note 
that  a  one-to-one  correspondence  between  prototype  and  model  temperature 
differences  in  each  equation  will  exist  only  when 


h*L*^ 

q* 


k  *L* 
w 

q* 


1 


(19) 


6 


from  which 


h*L* 

k  * 
w 


-  Bl* 


1 


(19a) 


where:  Bi  is  the  Biot  number* 

Bi*  »  Bi  /Bi 
m  p 

Ttie  requirement  that  Bi*  »  1  severely  restricts  the  flexibility  of  the 
model  to  simulate  prototype  values  of  h  and  k  ;  therefore,  this  require^ 
menr  cannot  always  be  met  and  temperature  distortions  will  result. 

There  are  two  methods  of  accounting  for  these  distortions;  (1)  imposing 
a  one-to-one  correspondence  between  model  and  prototype  shelter  wall 
temperatures,  and  (2)  imposing  a  one-to-one  correspondence  in  heat  flux 
betwetP  the  model  and  prototype 

Case  1  Condition 

Considering  the  Case  1  condition,  it  is  noted  from  Equation  15  that 
K^*L*  ■  q*.  Substituting  this  requirement  into  Equation  18, 


(t  -t  ) 
a  w'p 


h*L*^ 

k  *L* 
w 


(t  -t  ) 
a  w  m 


(18a) 


provided  ,  =  t 

w  w 
o  0 
p  m 


Since  Bi  *  hL/k 


w 


(t  -t  )  =  Bi*(t  -t  ) 

a  w'p  a  w  m 


(18b) 


V/har.  Bi*  1,  an  air  temperature  distortion  between  prototype  and  model 
will  occur. 


*The  Biot  number  has  the  same  grouping  of  variables  as  the  Nusselt 
number.  The  fundamental  difference  between  these  two  numbers  is  that 
the  Biot  number  involves  the  thermal  conductivity  of  a  solid,  and  the 
Nusselt  number  the  thermal  conductivity  of  a  fluid.  Thus  the  Blot 
number  is  an  index  of  the  relative  resistance  to  heat  flow  in  a  solid 
and  fluid,  whereas  the  Nusselt  number  expresses  a  comparison  of  the 
convective  capability  to  the  conductive  capability  of  the  fluid  only. 
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I 


(20) 


Let  At  ■  t  ~  t 

a  a  a 


Since  t  *  t  (by  dtfinition  for  Case  1) 

V  V  ■' 

m  p 


At  -  (t  -t  )  -  (t  -t  ) 

a  a  w^'m  a  w^p 

(21) 

■  <hA 

(22) 

(18b) 

(23) 

Substituting  from  Equations  22  and  23  into  Equation  21: 


At  *  “  (-^) 
a  Bi* 


(-^) 

Wp 


Or,  the  air  temperature  distortion  may  be  expressed  as 


(21a) 

(21b) 


At 


-) 


'a  Wp'  Bi* 

(Note  that  the  distortion  equals  zero  when  Bi*  »  1) 


(21c) 


The  mode 1> pro to type  air  temperature  relation  is  found  fiom  Equations  20 
and  21c,  where: 


-  (--^)  ( 


1-Bi*. 
Bi*  ^ 


(24) 


Since  the  model  may  not  have  the  same  initial  temperature  as  the 
chosen  prototype,  all  air  temperatures  should  be  replaced  by  temperature 
differences  between  temperatures  at  time  9  and  initial  temperatures  at 
time  0*0. 


,1-Bi*. 


(t  -t  )  »  (t  -t  )  -  (^)  (^i^-i^) 

a^  a_  p  a^  a  'm  '•hA'p^  Bi*  ' 


0 


*0 


(24a) 
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Case  2  Condlti jn 


If  the  Case  2  condition  Is  Imposed  (heat  flux  in  model  and 
prototype  are  equal),  the  model-prototype  air  temperature  relation  is 
found  to  be 


(t  -t  ) 

*0  %  P 


(t 


t  ) 


pm  0  0 


(See  Appendix  B  for  derivation) 

Note  that  the  Case  2  condition  Involves  model  distortion  of  both  the 
air  and  wall  temperatures,  whereas  the  Case  1  condition  results  in  model 
distortion  of  the  air  only.  Also  note  that  Case  1  and  Case  2  are 
separate  experimental  conditions.  Thus  Equations  24a  and  B-15a  are  not 
strictly  interchangeable.  Although  both  Cases  1  and  2  were  attempted 
experimentally,  only  Case  1  will  be  developed  in  detail  within  this 
report. 


Description  of  the  Model 

The  underground  shelter  model  was  a  1/7  scale  steel  arch  structure 
similar  to  that  used  for  several  actual  shelter  tests  which  have  been 
conducted. 3 See  Figure  2.  The  model  was  buried  in  dry  beach  sand  to 
assure  a  homogenous  soil  medium  and  reasonably  consistent  degree  of 
compaction,  and  also  to  eliminate  the  problem  of  moisture  migration. 

The  moisture  content  of  soil  affects  its  bulk  thermal  properties; 
thus  moisture  migration  during  a  test  would  result  in  a  requirement  to 
account  for  varying  thermal  properties.  (Effects  of  moisture  migration 
would  be  considerably  more  significant  for  the  model  than  the  prototype 
becruse  the  unsealed  model  temperature  gradients  in  the  soil  are  steeper.) 
The  thermal  properties  of  the  sand  (conductivity  and  specific  heat)  were 
measured  by  a  specially  constructed  conductivity  tester  and  standard 
calorimetric  techniques  respectively.  (The  conductivity  tester  was 
designed  to  measure  thermal  conductivity  by  means  of  steady-state  conduc¬ 
tion  between  two  concentric  cylindrical  surfaces.)  Model  temperatures 
were  measured  by  copper-constantan  thermocouples  placed  inside  the  shelter, 
on  the  shelter  walls  and  floor,  and  in  the  sand  surrounding  the  shelter. 
Thermocouples  for  measuring  soil  temperatures  were  soldered  to  brass  rings 
which  were  mounted  on  wooden  dowels  protruding  radially  from  the  shelter 
walls  and  floor  (see  Figure  2).  Heating  pad  units  stripped  of  padding 
were  used  in  conjunction  with  a  wattmeter  and  variac  to  provide  controlled 
sensible  heat  to  the  model.  The  pads  were  operated  far  below  their  rated 
output  to  minimize  radiant  heat  exchange.  Latent  heat  input  (part  of 
human  heat  output)  was  not  provided  for  the  model;  nor  were  a  ventilation 
system  and  an  air-conditioning  system  provided.  Instead,  the  sensible 
heat  output  was  modified  in  accordance  with  a  heat  balance  which  is 
described  below.  Although  this  procedure  resulted  in  an  air  temperature 
distortion  (which  will  be  discussed),  it  greatly  simplified  the  model 
system  and  allowed  maximum  experimental  control. 
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Heat  Input  to  the  Model 

The  heat  balance,  which  determined  the  heat  Input  to  the  model, 
required  consideration  of  human  sensible  and  latent  heat  outputs, 
climatic  conditions  for  the  geographical  location  of  the  prototype 
shelter  In  question,  and  ventilation  rate  and  alr-conditloning  capacity 
for  the  prototype  operating  conditions  being  modeled.  Sensible  heat 
was  added  to  the  model  according  to  the  net  heat  load  which  would  be 
Imposed  on  the  prototype  shelter  walls.  The  capacitance  effect  of  the 
shelter  air  could  be  neglected  because  the  thermal  capacitance  of  air 
compared  to  that  of  the  shelter  walls  Is  very  small. 

The  sensible  heat  load  Imposed  on  the  prototype  shelter  wall  Is: 


,  -  (835-8t^)  -  R(F)  - 

p(per  person) 

Human  heat  Air-  Ventilation 

output*  conditioning 

The  latent  heat  load  Imposed  on  the  prototype  shelter  wall  Is: 

q  -  (8t  -A35)^  -  R(l-F)  -  P„QH(w  -w  ) 

p(per  person) 

Human  heat  Air-  Ventilation 

output*  conditioning 

The  total  heat  load  Imposed  on  the  model  shelter  wall  Is: 


(25) 


(26) 


P  P  P 

when  q  ^0 

L 

p 

(27) 

qT  -qg 

when  qj^  <0 

(28) 

Notice  that  q  must  be  positive  to  be  considered,  A  latent  heat  load  Is 
."P 

manifested  by  moisture  condensation  on  the  prototype  shelter  walla. 
Evaporation  from  the  shelter  wall  requires  the  presence  of  moisture.  This 
mass  transfer  phenomena  could  not  be  controlled  analytically  for  the  model, 
thus  a  blanket  restriction  that  evaporation  from  the  shelter  walls  not  be 
considered  was  necessary. 

The  relative  humidity  Inside  prototype  shelters  Is  usually  very  high, 
thus  it  can  be  assumed  that  the  shelter  air  is  at  or  near  saturation  for 
purposes  of  computing  q^  .  and  q^  can  be  evaluated  from  Equations  25 

P  P  P 

*Mathematlcal  expression  suggested  by  Drucker^  to  approximate  metabolic 
heat  output  of  average  persons  engaged  in  moderate  activity.  (See  ASHRAE 
Guide  and  Data  Book,  Fundamentals  and  Equipment,  1961,  pp.  106-107.) 
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and  26  as  a  function  of  dry-bulb  temperature.  The  basic  problem  is  to 
assure  that  the  heat  input  to  the  model  is  varied  properly  according  to 
the  dry-bulb  temperature  inside  the  model.  At  this  point  it  should  be 
recalled  that  the  latent  heat  input  is  added  as  sensible  heat  to  provide 
greater  control  over  the  model.  Since  sensible  heat  transfer  between  the 
shelter  air  and  walls  r«i!quires  a  temperature  difference,  qg  -  hA(t^-t^), 

the  addition  c  ^  latent  heat  by  means  of  sensible  heat  results  in  another 
air  temperature  distortion.  This  distortion  can  be  corrected  as  follows: 

Let  (t;-t^)  -  (j^)p  -  (  hA  )P 

(29) 

(30) 

(tn-t  )  -  (^) 
aw  hA  p 

(31) 

since  1$  “II  ■  '•t 

P  P  P 

(27) 

(32) 

ta  +  -  2t^  " 

(52a) 

t'  ■  t  +  ) 

a  a  aw 

(32b) 

Let  At'  -  (t''-t  ) 

a  ■ 

% 

Then  At'  ■  (rr)  (t*'*  temperature  distortion) 
a  hA  p 

(3la) 

And 

(34) 

Incorporating  this  correction  into  Equations  24a  and 

B-13a  we  have: 

('.,-.,)p  •  -  4)p  -  <S>p  ‘ 

WV)  (24b) 

- 

<vv-  ■ ' VV-  ■  Vv’  * 

a  0 

(B-I5b) 
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Preparations  for  a  Model  Run 

Model  runs  were  made  using  both  the  Case  1  and  Case  2  conditions. 
However,  the  Case  1  condition  (which  imposes  equal  prototype  and  model 
wall  temperature  rises)  is  considerably  less  complex;  thus  this  condition 
will  be  used  exclusively  to  describe  the  procedure  for  setting  up  a  model 
run. 


As  mentioned  earlier,  the  basic  problem  of  setting  up  a  model  run 
is  properly  relating  heat  input  to  the  air  temperatures  inside  the  shelter. 
Once  the  climatic  conditions  at  the  chosen  geographical  location  have  been 
determined  from  weather  data,  and  prototype  values  for  the  many  physical 
and  operating  conditions  have  been  assigned,  the  heat  balances  (Equations 
25  and  26)  are  sol-  ed  for  the  range  of  air  temperatures  expected  to  occur 
in  the  prototype  shelter.  The  family  of  solutions  for  q^  as  a  function 

P 

of  t  are  the  basic  heat  inputs  to  the  prototype  shelter  wall.  These 

heat  inputs  must  be  modified  for  the  model  as  follows:  Recalling  that 
(t^-t  )  ■  ^^w”^w  ^  Case  1  condition.  Equation  15  requires  that 

^*0  ^  o  *** 


q*  ■  k  *L* 
^  w 


Since  q*  ■  q  /q 
m  p 


q  ■  k  *L*q  on  a  per  person  basis, 
m  w  p 


The  total  heat  input  to  the  oKxiel  is 


<It  "  "V^***! 

m  p 


where;  n  is  the  number  of  shelter  occupants. 


(35) 


(36) 


Now,  the  modal  air  temperatures  will  be  distorted  because  Bi*  my  not 
equal  unity,  and  the  latent  heat  portion,  q  ,  of  the  total  heat  input, 

i« 

P 

is  supplied  to  the  model  as  sensible  heat.  The  distortions  for  the 
Case  1  condition  are  part  of  Equation  24b.  These  distortions  are  found 
using  the  solutions  of  the  heat  balances  so  that  for  each  prototype 
shelter  air  temperatures,  t  ,  there  is  a  distorted  model  shelter  air 

*P 

tMq>erature  rise,  (t  -t  )  ;  and  the  model  heat  input,  q  ,  is  plotted 

9  o  m 

as  a  function  of  model  air  temperature  rise  (see  Figure  3).  Thl^  curve 
is  used  during  the  model  run.  After  the  run  has  been  completed,  model 
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air  temperature  distortions  are  applied  (with  a  change  in  algebraic  sign) 
to  the  experimental  air  temperature  data  to  yield  corrected  air  tempera¬ 
tures  for  the  prototype.  For  this  purpose  it  is  convenient  to  plot 
prototype  air  temperatures  as  a  function  of  model  air  temperature  rise 
(see  Figure  4).  Using  the  Case  1  condition,  the  model  gives  the  prototype 
shelter  wall  temperature  rise  directly  —  no  correction  is  necessary. 


Since  experimental  data  is  taken  in  real  time,  the  time  of  each  model 
temperature  observation  must  be  converted  to  prototype  time  after  the  run 
has  been  completed  --  the  equivalent  number  of  prototype  hours  for  each 
hour  of  model  (or  real)  time  must  be  found.  This  conversion  is  made  from 
Equation  10. 


0 

P 


0 

m 


(10) 


/  2 

The  conversion  factor,  /  L*  ,  Is  found  before  the  model  run  is  begun 
so  that  the  duration  of  the  run  can  be  anticipated. 


A  Typical  Model  Run 


To  examine  the  behavior  of  the  model  let  us  consider  a  model  run 
for  Washington,  D.  C.,  at  summer  climatic  conditions.  The  performance 
of  the  model  can  be  compared  to  analog  computer  results  for  the  same 
location  Ly  using  the  input  data  reported  by  Drucker.^  The  basic  climatic 
and  shelter  operating  conditions  were  evaluated  as  follows: 


Average  air  temperature  outside  shelter: 
Average  sir  absolute  humidity  outside  shelter: 
Initial  soil  temperature: 

Ventilation  rate  per  person: 

Air-condltlontng  rate  per  person: 


t^  -  7:‘>F 
w^  -  0.0162  Ib/lb 
t  es'^F 

w 

o 

Q  -  6  cfm 
R  -  0 


All  other  input  data  and  detailed  calculations  for  setting  up  the  model 
run  are  given  in  Appendix  C.  The  heat  input  curve  is  given  in  Figure  3, 
and  the  temperature  conversion  curve  Is  given  in  Figure  4. 

During  the  model  run,  teaq>eratures  were  recorded  as  a  function  of 
real  time.  Model  shelter  air  and  wall  temperatures  weie  determined  by 
averaging  the  response  of  several  thermocouples  in  each  instance.  As  the 
air  tcsq|>erature  changed,  the  heat  input,  ,  was  varied  according  to  the 

m 
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curve  shown  in  Figure  3.  After  the  experimental  data  had  been  taken, 
prototype  air  temperatures  were  found  from  the  distorted  model  air 
temperatures  by  means  of  Figure  4.  As  mentioned  earlier,  it  was  not 
necessary  to  correct  the  shelter  wall  temperature  data  (other  than  to 
account  for  the  difference  between  the  initial  soil  temperatures  of  the 
prototype  and  model).  Real  time  was  converted  to  prototype  time  by  means 
of  the  time  conversion  factor,  a  *^1,*^.  This  conversion  factor  was 
13.6;  thus  1  hour  of  model  time  was  equivalent  to  13.6  hours  of  prototype 
time  (see  Appendix  C). 

Results  of  the  Model  Run 

The  results  of  the  model  run  for  Washington,  D.  C.  are  shown  in 
Figures  5  and  6.  Figure  5  shows  the  results  when  using  the  Case  1 
condition,  and  Figure  6  shows  the  results  for  the  Case  2  condition.  The 
solid  line  curve  for  the  air  temperatures  were  developed  from  the  analog 
computer  results  obtained  by  Drucker.^  Reference  2  only  shows  the  analog 
results  for  the  shelter  conditions  at  the  tenth  day,  however.  Professor 
Drucker  was  kind  enough  to  provide  the  complete  analog  transient  solution 
for  purposes  of  comparison.  The  solid  line  shelter  wall  temperature  curve 
was  obtained  from  the  convection  equation  using  values  of  q^  calculated 

P 

from  Equation  25.  It  should  be  noted  that  the  model  results  using  both 
Cases  1  and  2  agree  closely  with  the  analog  computer  results.  It  should 
also  be  noted  that  the  transient  response  closely  resembles  a  step 
function.  The  air  temperature  rise  is  very  steep  initially,  and  then  it 
levels  off  so  that  by  the  fourth  or  fifth  day  of  occupancy  the  response 
is  quite  flat,  although  still  rising.  This  step-like  response  was 
observed  in  all  model  runs,  but  it  was  less  pronounced  in  situations 
where  the  total  air  temperature  rise  exceeded  15^F. 

The  results  shown  in  Figures  3  and  6  are  the  average  shelter 
temperature  conditions.  Actually  a  periodic  air  temperature  variation 
inside  the  shelter  occurs  because  the  outside  air  temperature  fluctuates 
periodically  over  a  24-hour  interval.  The  maximum  daily  air  temperature 
is  found  by  adding  a  temperature  increment  to  the  data  shown.  This 
increment  is  a  furiction  of  the  daily  outside  air  temperature  range  and 
the  ventilation  rate.  Teiiq>erature  increments  for  various  daily  outside 
air  temperature  ranges  and  shelter  ventilation  rates  have  been  determined 
from  Drucker 's  data,  and  can  be  found  from  the  family  of  curves  presented 
in  Figure  7.* 


^he  curves  shown  in  Figure  7  are  valid  for  the  shelter  area/person  and 
soil  thermal  property  values  used  by  Drucker  (see  Appendix  C).  However, 
it  is  believed  that  the  curves  have  general  applicability  because  the 
effect  of  shelter  area  per  person  and  soil  thermal  properties  on  the  basic 
transient  response  is  much  more  significant  than  the  effect  of  these 
parameters  on  inside  air  fluctuations. 
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The  average  dally  temperature  range  In  Washington,  D.  C.  during  the 
suoiner  Is  (according  to  Weather  Bureau  data).  From  Flgu..e  7,  a 

ventilation  rate  of  6  cfm/person  under  these  conditions  will  result  In  a 
temperature  Increment  of  1.7°F.  Thus  the  maximum  shelter  air  temperature 
occurring  during  the  tenth  day  of  occupancy  Is  84.5  +  1.7  «  36.20F.  Now, 
dry-bulb  temperature  alone  Is  not  a  complete  Index  of  comfort;  but  dry- 
bulb  temperature  considered  with  respect  to  absolute  humidity  (or  wet-bulb 
temperature)  leads  to  an  Index  of  comfort,  called  effective  temperature. 

Effective  temperature  Is  an  empirically  determined  Index  of  the 
degree  of  warmth  perceived  by  the  human  body  on  exposure  to  different  com¬ 
binations  of  temperature,  humidity,  and  air  movement.  The  temperature  was 
determined  by  trained  subjects  who  compared  the  relative  warmth  of  various 
air  conditions  In  two  adjoining  rooms  by  passing  back  and  forth  from  one 
room  to  another.  The  results  of  these  tests  have  led  to  an  effective 
temperature  chart  (see  ASHRAE  Guide  and  Data  book,  1961,  p.  109).  For 
the  sake  of  convenience,  a  portion  of  this  chart  has  been  transposed  to  a 
modified  psychrometrlc  chart  (Figure  8).  The  effective  temperature  lines 
have  been  drawn  for  stlll-alr  conditions. 


To  determine  the  maximum  effective  temperature  which  would  occur  In 
the  shelter  being  studied,  the  maximum  absolute  humidity  must  be  found. 

For  this  purpose,  recall  Equation  25  --  the  heat  balance  used  to  calculate 
the  latent  heat  load  Imposed  on  v’zhe  shelter  walls. 


q  -  (8t  -435)^  -  R(l-F)  -  P,QH(w  -w  ) 

L  a  p  a  a  A  p 

P 


(26) 


If  q  Is  positive,  condensation  Is  occurring  on  the  shelter  walls  and 
L 


the  air  Is  saturated  at  the  wail  temperature.  If  negative,  q  Is 

^P 

considered  equal  to  zero  and  the  absolute  humidity,  w^  ,  can  be  found  by 
rearranging  ^uatlon  26.  p 


(8t  -435)  -  R(l-F) 

—  pJh 

P  • 


''A  ^ 


(38) 


To  find  w 


max 


evaluate  Equation  38  using  t^ 


lx 


max 


-  86.2  F 
P 

a  -  0 
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w  -  0.0162  +  - - T 

^max  2.84  x  10 

P 

-  0.0162  +  0.009 

-  0.0252  Ib/lb 

Turning  to  the  psy chrome trie  chart.  Figure  8,  it  can  be  seen  that 
t  "  86.20F,  anc*  w  ■  0.0252  Ib/lb  correspond  to  an  effective 

*max  *max 

P  P 

temperature  of  very  nearly  85®F.  The  results  reported  by  Drucker  for 
the  tenth  day  of  shelter  occupancy  are  t  ■  86.5°F  and  ET  •  85®F. 

max 

P 


Model  Studies  on  Shelter  Shape 

The  excellent  agreement  between  the  model  study  results  and  the 
analog  computer  solutions,  as  illustrated  by  the  run  for  Washington,  D.  C., 
suggests  that  the  shape  of  the  shelter  may  not  significantly  affect  its 
transient  response  (the  Drucker  study  assumed  one- dimens tonal  heat  flow). 

If  the  assumption  of  one-dimensional  heat  flow  is  sound,  the  shelter 
surfaces  can  be  laid  out  as  the  plane  face  of  a  semi-infinite  slab.  Under 
these  conditions,  the  diffusion  equation  (6)  can  be  simplified  to 

Ox  w 

For  the  semi- inf initc  slab,  Carslaw  and  Jaeger^  suggest  a  solution  of  the 
form: 

2 

t  ■  c  w  (40) 


%fhich,  upon  substitution  into  Equation  39,  leads  to  an  expression  for  '^he 
teaqMirature  dist.. Ibution  in  the  slab  (soil)  as  a  function  of  time  following 
a  step  change  in  heat  flux. 


fa  0  ^  ■ 


2 

'4a  9^ 

w 


2  c 


At  x  «  0,  Equation  41  yields  an  expression  for  wall  temperatures  as  a 
function  of  time. 
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(42) 


t  + 
w 

o 


Js, 

k  A 
w 


a  e 


% 


The  model  was  tested  against  the  one-dimensional  heat  flow  relations 
as  follows:  at  time  0  »  0“,  q  *  0,  and  the  model  was  in  thermal  equili¬ 
brium;  at  time  0  »  ()+■,  the  heat  flux  was  raised  to  a  constant  value  equal 
to  the  maximum  heat  flux  which  would  occur  in  a  prototype  shelter  with 
normal  occupancy.  During  the  test,  wall  and  soil  temper'' tures  were 
measured,  but  the  heat  rate  was  not  changed.  The  experimental  data  was 
then  compared  to  theoretical  curves  obtained  from  Equations  ^l  and  42 
evaluated  in  terms  of  the  model  (see  Appendix  C  for  model  parameter 
values).  Since  the  model  was  not  being  run  with  respect  to  a  prototype, 
therewere  no  temperature  distortions  or  time  conversions.  However,  the 
total  heat  input  to  the  soil  was  roughly  equivalent  to  that  experienced 
with  a  model  run  for  a  typical  prototype. 

The  theoretical  curves  and  experimental  data  for  the  transient 
response  at  the  shelter  wall  are  shown  in  Figure  9;  and  the  theoretical 
curve  and  experimental  data  for  the  teraperature  distribution  in  the  soil 
after  8  hours  are  given  in  Figure  10.  The  experimental  results  agree 
very  satisfactorily  with  the  theoretical  curves.  Thl'.  evidence  indicates 
that  shelter  shape  effects  are  probably  not  serious.  Although  polydimen- 
sional  diffusion  from  the  shelter  edges  and  arch  was  expected  (this  was 
seen  by  plotting  isotherms  in  the  surrounding  soil),  the  exponential 
nature  of  the  temperature  gradient  (Figure  10)  shows  that  the  major 
portion  of  the  thermal  energy  absorbed  by  the  soil  occurred  close  to  the 
shelter  wall  thereby  lessening  the  effects  of  polydimensional  diffusion. 
For  this  reason,  it  is  believed  that  the  polydimensional  diffusion  effect 
would  have  been  minor  even  if  the  shelter  had  been  a  parallclc-piped. 

Use  of  Model  Study  Relationships  for  Prototype  S.helier  Runs 

It  is  conceivable  that  the  underground  shelter  transient  response 
could  be  studied  using  a  prototype  shelter  and  model  study  relationships. 
For  example,  having  run  a  prototype  shelter  bvried  in  a  given  soil,  it 
may  be  desirable  to  find  what  the  transient  r'isponse  wuuld  have  been  if 
the  shelter  had  been  burled  in  a  different  sell.  Assuming  that  the  soil 
thermal  properties  can  be  represented  by  single  values,  parameters 


pw*,  and  c  *  can  be  found, 
w 

Other  paratreters  are: 

q'*  •  1 

where:  ( 

)*  -  (  )j/(  ) 

2 

h^  -  1 

Subscrip' 

1  refers  to  the 
prototype. 

actual 

L*  *  1 

Subscript 

2  refers  to  the 

prototype  with 

modified  soil  properties. 
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Since  the  ratio  of  heat  fluxes  is  equal  to  unity,  a  Case  2  condition 

exists,  and  Equation  B-15b  is  applicable.  This  equation  can  be  greatly 

simplified  because  h*  =  1,  L*  =  1,  and  latent  heat  is  not  added  as 

sensible  heat;  however,  a  wall  temperature  distortion  occurs  because 

k  *  1. 

w 


(t  -t  )  =  (t  -t  )  +  (k  *-l)  (t  -t  ) 

a„  a  2  '  a  '1  w  w_  w  1 

0  0  0  0  8  o 


(B-15c) 


Since  the  initial  air  temperature,  t  ,  is  synonomous  for  both  (1)  and 

61 

O 

(2),  Equation  B-15c  can  be  further  simplified  to 


=  t  (k  *-l)  (t  -t  ) 

^0,  ”  '^8  '^o 


(B-15d) 


For  k  *  *'1,  the  modified  soil  thermal  conductivity  is  less  than  that  for 
the  actual  prototype,  and  one  would  expect  that  higher  air  temperatures 
would  occur  in  the  shelter  if  it  were  buried  in  such  a  soil  (assuming 
D^*  =  1,  c^*  =  1),  Checking  Equation  B-15d,  it  is  seen  that,  for 

k^*  »1,  a  temperature  increment  is  added  to  the  actual  prototype  air 

temperature  data  to  obtain  air  temperatures  for  the  modified  prototype 

soil.  The  converse  is  true  for  k  *  «=!. 

w 

The  time  conversion  factor  for  the  prototype  buried  in  soil  having 
modified  thermal  properties  is  found  from  Equation  10,  which  can  be 
simplified  because  L*  *  1 . 


e 


2 


a  *9, 

w  1 


(10a) 


where:  a  *  *  (k/pc)  * 

w  'w 

It  should  be  recalled  that  for  the  model-prototype  case  the  time  conversion 
factor  involved  which  allowed  many  hours  of  prototype  time  to  be  run 
per  hour  of  model  (or  real)  time.  For  the  prototype-prototype  case,  the 
time  conversion  factor  can  be  less  than  or  greater  than  unity.  Practically 
speaking  the  prototype-prototype  time  conversion  factor  will  not  vary  far 
from  unity  --  rarely  exceeding  the  range  from  1/2  to  2.  Moreover,  this 
conversion  factor  will  have  little  effect  on  the  shelter  transient 
response  because  the  response  closely  resembles  a  step  function  (see 
Figures  5  and  6). 

Since  the  prototype  shelter  heat  input  to  the  wall  varies  naturally 
with  climatic  conditions  and  conditions  of  the  air  inside  the  shelter, 
there  is  no  need  to  consider  the  heat  balances  (Equations  23  and  26)  and 
the  problems  associated  with  heat  control.  On  the  other  hand,  the 
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functional  character  of  the  heat  input  appears  to  preclude  a  non- 
dimensional  treatment  of  climatic  conditions;  thus  a  given  prototype 
installation  cannot  be  testec  for  a  range  of  climatic  conditions 
greater  than  that  which  occurs  naturally  unless  the  temperature  and 
humidity  of  the  ventilation  air  is  controlled,  or  controlled  heat  inputs 
are  employed  as  in  the  model. 


CONCLUSIONS  TO  PART  I 

As  a  tool  for  obtaining  engineering  data,  the  model  shelter  is  not 
completely  satisfactory  because  considerable  effort  is  required  to  set 
up,  conduct,  and  analyze  each  run.  However,  the  model  study  proved  to 
be  of  considerable  value  in  terms  of  developing  a  method  for  the  analysis 
of  complex  heat  transfer  systems  and  generally  substantiating  the  computer 
programs  which  have  been  developed  by  others.  More  specifically  the 
results  of  model  experimentation  have  shown  that; 

1.  The  thermal  transient  response  of  an  underground  shelter  can  be 
examined  by  model  studies. 

2.  The  analog  computer  studies  conducted  by  Drucker,  et  al,^  for 
the  fully  buried  underground  steel  arch  shelter  are  probably  sound. 

3.  Shelter  shape  does  not  seriously  affect  the  thermal  transient 
response.  Shape  effects  include  temperature  distortions  due  to  poly- 
dimensional  heat  diffusion  at  the  edges  and  over  the  arch.  These  effects 
are  minor  because  the  temperature  gradient  does  not  extend  appreciably 
into  the  soil  surrounding  the  shelter  during  the  expected  period  of 
occupancy . 


4.  One-dimensional  conduction  in  the  soil  can  be  assumed  with 
confidence. 
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PART  II.  NON- COMPUTER  DESIGN  PROCEDURE  FOR  DETERMINING  THE  CAPACITY 
OF  UNDERGROUND  SHELTER  ENVIRONMENT  CONTROL  SYSTEMS 


Introduction 

The  non-computer  design  procedure  for  determining  the  capacity  of 
underground  shelter  environment  control  systems  depends  on  analytical 
tools  developed  in  Part  I.  Briefly  the  procedure  is  as  follows:  The 
shelter  soil  system  Is  conceived  as  a  semi- inf inite  slab  having  a 
plane  face  equal  in  area  to  the  total  shelter  area  (per  person  basis). 

The  maximum  dry-bulb  temperature  and  effective  temperature  tolerable 
(or  desired)  inside  the  shelter  at  the  end  of  the  expected  occupancy 
period  are  determined,  and  trial  values  for  the  ventilation  rate  and 
air-conditioning  capacity  are  assigned.  Next,  the  maximum  dry-bulb 
temperature  is  converted  to  an  average  value  to  account  for  daily  air 
temperature  fluctuations  inside  the  shelter.  It  is  assumed  that  the 
air  temperature  rises  to  the  average  value  upon  occupancy  and  remains 
at  this  value  during  the  entire  occupancy  period.  From  the  heat 
balance  equations  (25  and  26)  evaluated  for  the  average  Inside  air 
temperature  and  climatic  conditions  at  the  intended  shelter  location, 
the  sensible  and  latent  heat  rate  inputs  to  the  shelter  wall  are  deter¬ 
mined.  Now,  responding  to  the  rise  in  air  temperature,  the  wall 
temperature  rises  with  time  in  accordance  with  the  expression  developed 
for  a  semi- inf inite  slab  subject  to  a  step  Increase  in  heat  flux 
(Equation  42).  Eventually  a  maximum  wall  temperature  is  reached  above 
which  the  sensible  heat  to  the  shelter  wall  cannot  be  supported  at  the 
average  air  temperature,  and  the  occupancy  period  must  be  terminated. 

Once  the  maximum  wall  temperature  has  been  ascertained  (from  the  con¬ 
vection  equation),  the  allowable  occupancy  period  Is  determined  from 
Equation  42.  If  the  allowable  occupancy  period  Is  less  than  that 
required,  the  ventilation  rate  and/or  alr-conditloning  rate  must  be 
adjusted  to  lower  the  heat  Input  to  the  shelter  walls.  If  the  allowable 
occupancy  period  Is  considerably  greater  than  that  required,  the 
environmental  control  system  capacity  should  be  reduced. 

It  should  be  noted  that  two  fundamental  assumptions  are  necessary 
for  the  design  procedure:  (1)  the  air  temperature  response  In  an  actual 
shelter  can  be  thought  of  as  a  step  function,  and  (2)  one -dimensional 
heat  flow  from  the  walls  and  floor  of  an  actual  shelter  Into  the  surround¬ 
ing  soil  prevails.  The  bases  for  these  assumptions  are  covered  In  Fart  I. 
Other  assumptions  are  given  In  Appendix  A. 

The  design  procedure  requires  a  certain  degree  of  trial  and  error; 
however.  It  Is  straightforward,  flexible,  and  can  be  used  without  difficulty 
by  trained  personnel.  The  use  of  the  design  techniques  will  become 
apparent  from  the  step-by-step  procedural  outline  and  numerical  example 
below. 
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Outline  of  the  Design  Procedure 

I.  Determine  size  of  the  shelter  and  number  of  occupants  (usually 
governed  by  operational  requirements).  Ascertain  24-hour  average  cli¬ 
matic  conditions  (dry-bulb  temperature  and  wet-bulb  temperature)  and 

the  daily  temperature  range  at  the  intended  location  from  weather  data.^>7 
Entering  the  psychrometric  chart  with  dry-bulb  and  wet-bulb  temperature 
(or  relative  humidity)  data,  the  absolute  humidity  can  be  found. 

Determine  (or  assume)  soil  thermal  properties  at  the  intended  location. 

II.  Determine  maximum  dry-bulb  temperature  and  effective 
temperature  tolerable  (or  desired)  inside  the  shelter  at  the  end  of  the 
expected  occupancy  period  (usually  10  -  14  days).  The  following  table  of 
acceptable  and  tolerable  thermal  limits  for  healthy  persons  at  rest  and 
properly  clothed  is  suggested  by  Yaglou:® 


Effective  Temperature 
(minimum  air  movement) 

Lowest  acceptable  for  continuous  exposure- 
manual  dexterity  may  be  affected. 

50°F 

Optimum  for  comfort,  with  60%  relative 
humidity. 

68®F  -  72°F 

Perspiration  threshold.  Acceptable  for 
continuous  exposure. 

78°F 

Endurable  in  emergencies  for  at  least  two 
weeks—  possible  heat  rash  in  prolonged 
exposures. 

850F 

Possible  heat  exhaustion  in  unacclimatized 
persons. 

88°F 

Possible  heat  exhaustion  in  acclimatized 
persons . 

920f 

III.  Assuma  a  trial  ventilation  rate  and  air-conditioning  capacity. 

IV.  Test  the  ventilation  rate  and  air-conditioning  capacity  for  the 

maximum  air  temperature  by  ascertaining  effective  temperature  through 
consideration  of  the  latent  heat  balance  (evaluate  Equation  26  )• 

*nax 
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A.  Latent  heat  output  per  person. 


q,  -  (8t^  -435)  (43) 

max 

B.  Latent  heat  removed  by  air-conditioning,  per  person. 

q  -  (l-F)R  (44) 

^r 

C.  Latent  heat  removed  by  ventilation  air,  per  person. 

V  sat 

where:  w  is  the  absolute  humidity  of  air  saturated  at  the 

*8at 

maximum  shelter  air  dry-bulb  temperature,  t 

*max 


D.  Net  latent  heat  to  the  shelter  wall. 


(46) 


A  positive  value  of  q.  indicates  a  saturated  air  condition  at  t 

w  max 

and  the  effective  temperature  will  equal  t  If  this  effective 

max 

tMsperature  exceeds  the  effective  temperature  allo%red  in  Step  11,  the 
ventilation  rate  and/or  air-conditioning  capacity  should  be  increased. 
If  q  is  negative,  it  is  considered  equal  to  aero  (see  Section  1,  Heat 

V 

Input  to  the  Hodel),  and  w^  must  be  calculated  from  Equation  38a. 


*'.  *  "a  +  -pTpr^ 

Now  enter  the  modified  psychrometrlc  chart  with  t  and  w  to  find  the 

max  * 

effective  temperature.  If  this  value  exceeds  the  maxlsnaB  effective 
temperature  allowed  in  Step  11,  the  ventilation  rate  and/or  alr-conditloninf 
capacity  should  be  Increased. 

When  the  envlronsMot  control  system  capacity  is  adequately  slsed  for 
t  and  ET^^,  the  allowable  duration  of  shelter  occupancy  under  these 

*max 

conditions  can  be  found  from  the  steps  following. 
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V.  Having  completed  Step  IV  so  as  to  satisfy  the  effective 
temperature  requirements  determined  in  Step  11,  find  the  average  dry- 
bulb  temperature,  t  ,  tolerable  (or  desired)  at  the  end  of  the  expected 
shelter  occupancy  period  by  subtracting  a  trap£.’ratu^‘e  increment  from 

t  This  temperature  increment,  which  accounts  for  daily  temperature 

*max 

fluctuations  inside  the  shelter,  is  obtained  from  Figure  7.  Note  that 
the  increment  is  a  function  of  ventilation  rate  and  daily  temperature 
range  at  the  intended  location. 

VI.  Now  calculate  the  net  heat  input  to  the  shelter  walls  by 
evaluating  the  heat  balance  equations  using  t^. 

A.  Sensible  heat  balance 

1.  Heat  generated  per  persot? 

Og  -  (ft35.8t^)  (47) 

g 


2.  Heat  removed  hy  air-conditioning,  per  person 

qg  -RF 


(A8) 


3.  Heat  removed  by  ventilation  air,  per  person 


«ihere: 


V 

4.  Net  sensible  heat  to  the  shelter  wall. 

^s 

w  g  r  V 

B.  Latent  heat  balance  (follows  format  of  Step  IV) 

q,  -  (8t^-  435)  -  (l-F)R  -  P^QH(w^  -  w^) 
w  sat 

w  is  the  absolute  humidity  of  air  saturated  at  t  , 

*sat  * 

C.  Total  net  heat  input  to  the  shelter  wall. 


•t  ’’s 

WWW 


‘‘t  ’“s 


wham 


when 


q,  so 

S# 


(49) 


(50) 


(26) 


(27) 

(28) 
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From  Section  1,  q.  is  restricted  to  positive  values  because  evaporation 

L 

V 

from  the  shelter  walls  is  difficult  to  manage  analytically. 

Vll.  Having  determined  the  heat  inputs  to  the  shelter  wall  for  a 
given  air  temperature,  t  ,  find  the  wall  temperature  required  to  satisfy 
the  convection  equation  for  sensible  heat  transfer.  This  wall  tempera¬ 
ture  is  the  maximum  value  allowed  at  the  end  of  shelter  occupancy. 


t 

w 

max 


(51) 


where:  A  is  the  total  shelter  area  per  person. 

Vlll.  Finally,  calculate  the  allowable  duration  of  shelter 
occupancy  under  the  set  conditions.  The  tra.^sient  response  of  a  semi- 
infinite  slab  subject  to  a  step  increase  in  heat  flux  has  been  considered 
in  Part  1,  Model  Studies  on  Shelter  Shape.  At  the  wall,  the  temperature 
response  is  expressed  by  Equation  42.  This  equation  can  be  rearranged 
to  give  the  allowable  duration  of  shelter  occupancy  (when  t  is  reached). 

max 


Simplifying 


2 


hours 


e 


^•9Cq^ 

*w 


«fhere:  t  is  the  initial  soil  temperature 
o 


days 


(52) 


(53) 


(54) 


C 


AV(k^ 
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If  the  result  of  Equation  (54)  is  equal  to  or  greater  than  the  expected 
duration  of  shelter  occupancy,  the  ventilation  rate  and  alr-conditloning 
capacity  assumed  in  Step  III  are  sufficient.  If  the  calculated  duration, 
d,  greatly  exceeds  the  expected  occupancy  term,  the  ventilation  rate 
and/or  the  alr>condltlonlng  rate  should  be  reduced  to  provide  a  more 
economical  environment  control  system.  Naturally,  If  the  calculated 
duration  Is  less  than  the  expected  occupancy  period,  the  capacity  of  the 
environment  control  system  should  be  Increased. 

Numerical  Example 

As  a  numerical  example  of  the  hand  calculation  method,  consider  the 
shelter  and  conditions  used  to  Illustrate  the  model  study.  Let 
t^  «  86°F  and  ET  *  85°F  at  the  end  of  shelter  occupancy  (results 
max 

of  the  model  study).  Recall  that  Q  «  6  cfm/person  and  R  ■  0.  Further* 
more,  t.  ■  77°F,  t  ■  68®F,  w  *  0.0162  Ib/lb,  and  the  dally  temperature 

A  W 

O 

range  Is  18^.  See  Appendix  C  for  the  soil  properties  and  other  parameter 
values.  Having  completed  Steps  I  •  III,  the  trial  ventilation  rate  and 
alr-conditloning  capacity  are  tested. 


Step  IV.  Determine  ET  from  latent  heat  balance  evaluated 

max 


q.  -  (8x86-435)  -  2.84  x  10^(0.0272-0.0162) 
w 

•  -59  Btu/hr-person 

Since  q.  Is  negative.  It  Is  considered  equal  to  zero,  and  w  must 

L  ft 

V  Q&X 

be  found  from  Equation  38a. 


w  .  0.0162  4 
•max  2.84x10 

-  0.0251  Ib/lb 

Entering  the  modified  psychrometrlc  chart  (Figure  8)  with 
t  ■  86°r  and  w  ■  0.0251  Ib/lb,  ET  is  found  to  be  Just  slightly 

ft  ^  ft  flUUC 

mftx  «ftx 

less  than  85^F;  thus  the  ventilation  rate  and  alr-conditloning  capacity 
ap^ar  to  be  adequate  for  t^  However,  the  duration  of  allowable 

max 

shelter  occupancy  for  these  conditions  must  still  be  found. 
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Step  V.  The  average  shelter  air  temperature,  t  ,  at  the  end  of 

tt 

10  days  is  found  by  subtracting  an  increment,  At  ,  from  t  .  Entering 

max 

Figure  7  with  a  daily  temperature  range  of  18°F  and  ventilation  rate  at 
6  cfm/person.  At  *  1.7°F,  and  thus  t  «  (86-1.7)  *  84.3°F. 

A  A 

Step  VI.  Net  heat  input  to  shelter  walls. 

A.  Sensible  heat 

qg  =  (835-8x84.3)  -  6.5(84.3-77) 
w 

■  113.6  Btu/hr-person 


B.  Latent  heat 

q,  -  (8x84.3-435)  -  2.84x10^(0.0258-0.0162) 

Ij 

w 


■  -33  Btu/hr-person 

C.  Total  heat 


since  q 


w  w  w 

■  113.6  Btu/hr-person 


Step  VII.  Maximum  wall  temperature 


w 


max 


-  84.3  - 


-  81.6  F 


113,6 

42.4 


;o. 


Step  VIII.  Allowable  duration  of  shelter  occupancy. 


9 


(81.6-68) 

73 - - 

4.9x7,85x10  xll3.6 


■  10.1  days 


Since  the  expected  shelter  occupancy  period  for  this  case  is  10  days,  the 
ventilation  and  air-conditioning  rates  appear  to  be  sufficient.  It  should 
be  noted  that  the  non-computer  design  procedure  has  given  results  vrhlch 
agree  very  closely  with  the  model  study  and  analog  computer  results 
reported  in  Part  I. 
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Discussion 


When  the  hand  calculation  method  is  used  for  a  number  of  assumed 
air  temperature  (t  )  conditions,  the  family  of  solutions  generated  takes 

on  the  appearance  of  a  tranalet t  response  even  though  the  solutions  do 
not  actually  constitute  a  transient.  The  family  of  solutions  for  the 
particular  shelter  considered  in  this  report  is  shown  in  Figure  11.  It 
should  be  noted  that  significant  temperature  discrepancies  between  the 
hand  calculation  solutions  and  the  Drucker  analog  solution  (true  transient) 
occur  at  duration  times  less  than  five  days,  but  that  the  two  solutions 
appear  to  converge  at  approximately  10  days.  Since  we  are  really  only 
interested  in  the  temperature  condition  of  the  shelter  at  the  end  of  the 
expected  occupancy  period,  the  aforementioned  temperature  discrepan^'les 
are  of  small  consequence. 


As  a  matter  of  interest,  let  us  consider  the  nature  of  the  hand 
calculation  method  to  appreciate  the  cause  of  the  temperature  discre¬ 
pancies  seen  in  Figure  10.  To  begin  with,  the  hand  calculation  method 
is  similar  to  the  model  study  analysis  in  that  the  net  heat  transfer  to 
the  shelter  walls  is  found  by  a  heat  balance.  However,  when  using  the 
hand  calculation  method,  the  shelter  air  temperature  response  is  idealized 
to  a  step  function  when  the  final  air  temperature,  t  ,  is  assumed  in 

*max 

Step  II.  (Remember  that  this  temperature  is  reduced  to  an  average 
temperature  which  is  held  constant  over  the  entire  shelter  occupancy 
period.)  The  use  of  a  constant  air  temperature  is  required  to  fix  a 
constant  heat  input,  q^  ,  to  the  shelter  walls,  and  thereby  greatly 

w 

simplify  the  calculation  procedure.  If  q^  were  allowed  to  vary  with 


time,  a  finite  difference  solution  would  be  necessary;  such  a  solution  is 
more  appropriate  for  a  computer  than  hand  calculation.  Once  the  constant 
heat  flow  to  the  shelter  wall  has  been  established,  the  wall  temperature 
rises  according  to  Equation  42  to  a  temperature,  t  ,  determined  by 


Equation  51  which  was  derived  from  the  convection  equation.  At  this 
point,  the  equations  governing  the  heat  transfer  phenomena  are  satisfied 
and  a  solution  exists.  The  time  at  which  a  solution  exists  for  a  given 
heat  transfer  rate  to  the  shelter  wall  Is  the  upper  bound  for  shelter 
occupancy  (duration)  at  the  assumed  t  Now,  the  hand  calculation 

*max 

method  works  because  a  prototype  shelter  exhibits  an  extresMly  sharp 
initial  air  tmsperature  rise  which  approximates  a  step  function  response. 
However,  remember  that  heat  flow  to  the  shelter  wall  Is  a  function  of 
shelter  air  temperature.  Thus,  the  extent  to  vhlch  a  prototype  shelter 
deviates  from  the  Idealized  step  function  response  will  determine  the 
degree  of  discrepancy  between  the  hand  calculation  solution  and  prototype 
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t 


behavior.  Kt  tlnte  days,  the  prototype  temperature  response  over 

the  total  time  Interval  shows  a  significant  deviation  from  a  step 
function  respon5e;  but  at  time  G  »  10  days,  the  temperature  response  is 
flat,  and  over  the  tot?l  time  inter^'a!  the  initial  deviation  from  a  step 
response  is  less  signif icai-.t  (see  Figure  11). 

The  ha. id  calculation  solution  can  be  applied  to  the  buttoned-up 
shelter  situation  without  theoretical  temperature  discrepancies  from 
prototype  behavior  because  heat  flow  to  the  shelter  wall  is  not  a 
function  of  shelter  air  temperatiire.  The  "ventilation"  terms  of 
Equations  46  and  30  which  cause  the  functional  relationship  between 
air  temperature  and  heat  flow  are  no  longer  considered.  In  this  case; 

It  •  ’s  ’l 

w  w  w 

»  (835-8t  )  -  Rf  +  (8t  -  435)  -  R(l-F)  (66) 

A  A 

*  400  -  R  Btu/hr-person 

When  there  is  no  air-conditioning,  the  heat  flow  to  the  shelter  wall  is 
simply  equal  to  the  total  metabolic  heat  output  of  the  occupants. 

Additional  Considerations  for  the  Design  Procedure 

Although  the  envirojimeut  control  system  for  the  shelter  located  in 
Washington,  D.  C.  appears  to  be  sufficient  fo;  the  chosen  design  condi¬ 
tions,  the  conditions  should  be  exsnined  to  determine  limitetions  of 
shelter  habltebility . 

First  of  all,  the  chosen  everege  outdoor  temperature  of  77°F  for 
Washington,  D.  C.  represents  a  lOX  design  condition.  A  temperature 
expressed  In  terms  of  X  design  condition  represents  s  climetic  condition 
which  prevails  or  is  exceeded  during  the  summer  months  (June  through 
September).  For  example,  the  lOX  design  dry-bulb  everege  temperature 
corresponds  roughly  with  the  24-hour  average  cemperetures  during  the 
suBttcr  ssonths.  Raising  the  lOX  design  dry-bulb  temperature  for  a  given 
location  by  5°F  usually  leads  to  design  conditions  in  the  I  to  S*  range, 
depending  on  location.  For  Washington,  D,  C.,  the  4%  design  dry-bulb 
teaq>ereture  Is  62<^F,  end  the  IX  design  dry-bulb  temperature  is  86^F. 

The  point  le  that  the  24-hour  everege  tempereture  In  Washington,  D.  C. 
could  exceed  77^F  for  several  consecutive  days,  end  a  shelter  having  e 
merglnel  environmental  control  system  designed  for  lOX  conditions  could 
become  dangerously  uncomfortable.  This  argument  also  extends  to  the  use 
of  wet-bulb  temperature  design  date. 

Another  point  which  deserves  consideration  is  the  Mvebolic  heat 
output  expressions  used  in  this  study  (part  of  Equations  23  and  26) 
which  are  based  on  s  total  heat  out^>uc  cf  400  8tu/hr-person.  Total 
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metabolic  heat  c-  'cput  depends  o^  the  general  activity  level  of  shelter 
personnel.  Although  the  total  value  of  AOO  Btu/hr-person  is  widely 
accepted,  Strope,  et  al,^  recoromend  an  average  total  heat  output  of 
approximately  500  Btu/hr-person.  If  this  recommended  value  had  been 
used  for  the  shelters  considej.ed  in  this  report,  the  required  capacity  of 
the  environment  control  system  would  have  been  somewhat  greater.  The 
determination  of  sensible  and  latent  heat  outputs  as  a  function  of  air 
temperature  for  total  heat  c  itputs  exceeding  400  Btu/hr-person  can  be 
found  from  the  ASHRAE  Guide  and  Data  Book.^ 

Additional  Solutions  Using  the  Design  Procedure 

Following  the  work  of  Drucker,^  the  non-computer  design  procedure 
was  applied  to  various  shelter  conditions  at  a  number  of  locations 
where  summer  climatic  conditions  might  be  severe  enough  to  cause  shelter 
habitability  problems.  Tne  solutions  are  summarized  in  tabular  form 
below.  These  solutions  should  not  be  used  for  design  purposes  because 
the  non-computer  design  'procedure  has  not  been  fully  substantiated; 
however,  they  are  useful  for  comparison  purposes  to  determine  the  effect 
of  climate,  shelter  population  level,  and  Initial  soil  temperature.  The 
shelter  physical  parameters  and  soil  properties  remain  unchanged  from 
those  values  given  in  Appendix  C.  Maximum  effective  temperature  allowed 
was  850F. 
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ST.  LOUIS,  MISSOURI 
(Mississippi  Basin) 
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Conclusions 


1.  Based  on  a  review  of  the  non-computer  design  procedure,  the  most 
Influential  parameters  affecting  the  capacity  requirements  for  environ¬ 
ment  control  systems  appear  to  be: 

a.  Total  metabolic  heat  output  per  person  and  proportion  of 
sensible  to  latent  heat. 

b.  Outdoor  climatic  conditions  (Including  daily  dry-bulb  temperature 
variations) . 

c.  Initial  soil  temperature. 

d.  Soil  thermal  conductivity. 

e.  Total  shelter  surface  area  per  person. 

2.  Comparing  the  results  of  the  additional  solutions  In  the  previous 
section,  a  number  of  trends  become  evident:  (valid  only  for  the  locations 
considered). 

a.  Ventilation  alone  at  3  cfm/person,  considered  adequate  for 
control  of  exygen  and  carbon-dioxide  levels,  Is  not  sufficient 
to  maintain  effective  temperatures  below  during  the  assumed 
shelter  occupancy  period. 

b.  For  the  10%  design  condition,  ventilation  alone  In  excess  of 
3  cfm/person  Is  sufficient  for  the  locations  studied  except 
Houston  and  Phoenix.  For  these  two  locations,  note  that  the 
lnlt:^al  soil  temperatures  and  average  air  temperatures  are 
relatively  high. 

c.  An  Increase  of  5^  In  the  outside  air  temperatures  (over  the 
dally  range)  from  the  10%  design  conditions  will  require  the 
use  of  air  conditioning  In  Washington,  Houston,  and  Phoenix. 
However,  In  Providence,  Chicago,  and  St.  Louis,  where  initial 
soil  temperatures  are  relatively  low,  no  air-conditioning  is 
required  If  ventilation  rates  are  Increased  over  that  required 
for  10%  conditions.  This  latter  situation  occurs  because  the 
soil  acts  as  a  heat  sink,  and  also  because  ventilation  air  is 
capable  of  carrying  away  considerable  amounts  of  latent  heat 
In  the  form  of  water  vapor. 
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d.  Under  17.  design  conditions,  the  use  of  air  conditioning  is 
required  in  all  locations  except  Providence  where  an  increase 
in  ventilation  rate  over  the  57»  condition  is  sufficient.  In 
Chicago,  St.  Louis,  and  Houston,  the  use  of  air  conditioning 
is  accompanied  by  a  reduction  in  ventilation  rate  from  the 

4  to  57.  design  condition.  This  reduction  in  ventilation  rate 
occurs  because  high  outdoor  dry-bulb  temperatures  (greater 
than  t^  allowed)  result  in  a  sensible  heat  gain  to  the  shelter. 

In  Phoenix,  the  ventilation  adds  sensible  heat  to  the  shelter 
even  at  the  10%  conditions,  thus  ventilation  is  maintained  at 
low  levels  for  all  design  conditions. 

e.  For  all  locations  except  Houston  and  Phoenix,  increasing 
shelter  population  by  30%  can  be  taken  care  of  under  the  10% 
design  weather  conditions,  by  increasing  the  ventilation  rate 
per  person  over  that  required  for  normal  occupancy;  the  venti¬ 
lation  rate  at  Houston  and  Phoenix  is  more  than  1507.  that  which 
is  normally  required.  Due  to  the  high  outdoor  dry-bulb  tempera¬ 
ture  at  Houston  and  Phoenix  the  ventilation  rate  must  be  held  to 
low  levels  and  the  air  conditioning  capacity  be  increased 
accordingly.  Shelter  overpopulation  has  the  effect  of  reducing 
the  shelter  surface  area  per  person. 

f.  If  a  24-hour  buttoned-up  period  is  required,  air  conditioning  is 
necessary  for  those  locations  where  the  initial  soil  temperature 
is  above  68°F  (Washington,  St.  Louis,  Houston,  Phoenix). 

These  conclusions  are  not  intended  as  design  rules;  indeed,  they  are 
not  specific  enough,  and  they  are  based  on  calculations  for  a  particular 
set  of  shelter  and  soil  parameter  values  (Appendix  C) .  The  intent  is  to 
demonstrate  that  a  shelter  environs^nt  control  system  can,  and  should, 
be  tailored  to  the  specific  location  of  the  shelter.  Continuing 
environment  control  studies  are  devoted  to  giving  the  shelter  designer 
effective  analytical  methods  to  meet  this  goal. 

Suggvtstions  for  Future  Work 

1.  Investigate  simplifying  the  non-computer  design  procedure  by  the  use 
of  charts,  graphs,  or  nomograms. 

2.  If  sufficient  data  from  an  actual  shelter  test  can  be  found,  compare 
the  results  of  the  test  with  a  solution  of  the  design  procedure  evaluated 
with  the  actual  test  parameter  values.  Two  mejor  problems  will  be  to 
determine  representative  soil  thermal  properties  and  24-hour  average 
climatic  conditions. 
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3.  Expand  the  non-computer  design  procedure  to  Include  concrete 
structures.  A  suggested  solution  method  is  to  use  an  equivalent 
coefficient  of  heat  penetration  (kpc)»  in  Equation  54. 

(kpc)  ,  »  A(kpc)  +  B(kpc)  . 

equiv.  '  v  r  concrete 

where:  A  and  B  are  coefficients  which  account  for  the  relative 

thermal  energy  absorption  of  concrete  and  soil  in  the 
composite  wall.  Both  A  and  B  are  positive,  and  A  +  B  *  1 

4.  Expand  the  non-computer  design  procedure  to  include  partially  burled 
structures.  In  this  case  only  part  of  the  shelter  surface  area  would  be 
considered  with  respect  to  the  semi-infinite  slab  (heat  sink).  The 
balance  of  the  area  would  be  considered  as  the  plane  face  of  a  thick 
wall  subject  to  a  24-hour  average  sol-air  temperature.  The  net  heat 
into  or  out  of  this  thick  wall  would  be  considered  in  determining  the 
heat  load  to  the  semi- infinite  slab. 

5.  Expand  the  design  procedure  to  Include  "basement”  shelters  in  which 
all  walls  contact  the  soil  (semi- infinite  slab).  In  this  case  the  net 
heat  load  through  the  celling  of  the  shelter  can  be  determined  by 
standard  air  conditioning  techniques  for  the  thermal  evaluation  of 
buildings.  This  net  heat  load  would  be  considered  in  determining  the 
heat  load  to  the  semi- Infinite  slab. 
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PART  III.  GENERAL  RECOMMENDATIONS 


1.  The  use  of  marginal  design  criteria  with  respect  to  controlling  the 
thermal  environment  of  shelters  should  be  carefully  studied.  An 
Inadequately  designed  shelter  may  offer  protection  from  fallout,  but 

It  may  also  offer  extreme  discomfort,  disease,  and  even  death  from 
Intolerably  high  thermal  levels. 

2.  As  a  step  toward  better  shelter  design,  the  metabolic  heat  output  of 
people  In  various  modes  of  activity  associated  with  normal  shelter  life 
should  be  better  understood. 

3.  The  design  of  the  shelter  thermal  environment  control  system  should 
be  based  on  a  careful  analysis  of  the  climatic  conditions  at  the  Intended 
location.  Whenever  possible,  conservative  design  data  (less  than  lOX) 
should  be  used. 

The  possibility  of  a  low-cost  packaged  air  conditioning  unit  suitable 
for  shelter  use  should  be  Investigated. 

5.  The  proper  management  of  shelter  ventilation  as  a  function  of  outside 
climatic  conditions  might  reduce  or  eliminate  the  air  conditioning  require¬ 
ments  In  some  locations  where  excessive  dry-bulb  temperatures  occur  during 
the  day. 

6.  An  economic  study  of  achieving  lowest  cost  environment  control  systems 
by  trade-offs  between  ventilation  and  air  conditioning  should  be  under¬ 
taken.  Considering  the  problem  of  collective  protection.  It  might  be 
more  economical  to  provide  air  conditioning  with  reduced  ventilation  rates. 
Moreover,  In  many  cases,  this  approach  will  adequately  control  the  shelter 
environment  up  to  1%  design  conditions. 
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Appendix  A 


ASSUMPTIONS 


".lie  Prototype  Shelter 

1.  The  shelter  is  sufficiently  buried  so  that  temperature  gradients  in 
the  soil  extending  from  the  surface  will  not  affect  the  temperature 
response  of  the  shelter. 

2.  The  initial  soil  temperature  is  uniform. 

3.  The  soil  is  homogenous,  and  its  thermal  properties  are  uniform. 

4.  The  soil  thermal  properties  are  independent  of  temperature. 

5.  The  shelter  metal  skin  and  concrete  deck  can  be  treated  as  part  of 
the  soil  system  without  introducing  serious  error. 

6.  The  thermal  properties  of  air  are  uniform  and  independent  of  temperature. 

7.  The  temperature  of  the  inlet  air  varies  sinusoidally  between  the 
daytime  maximum  and  nigh time  minimum. 

6.  The  absolute  humidity  of  the  inlet  air  is  constant. 

9.  Perfect  air  mixing  occurs  inside  the  shelter. 

10.  The  convective  heat  transfer  coefficient  Inside  the  shelter  is 
constant . 

11.  Sensible  and  latent  heat  generated  by  the  shelter  occupants  can 
expressed  by  the  relations:  (there  are  no  other  heat  sources^. 

q  •  (835-8t  )  Stu/hr-.^erson 

s  * 

q  ■  (8t  -435)  Btu/hr-person 
I#  ^ 

12.  Condensation  which  occurs  on  shelter  walls  docs  not  evaporate  even 
if  the  surface  temperature  Subsequently  rises. 

13.  The  effect  of  solar  radiation  upon  the  surface  of  the  soil  cen  be 
neglected . 


41 


I 


The  Model  Shelter 


1.  All  assumptions  for  the  prototype  shelter  extend  to  the  model 
shelter. 

2.  The  heat  response  of  the  heating  pads  is  instantaneous. 

3.  Net  radiation  from  the  heating  pads  to  the  she  . ter  wall  can  be 
neglected. 

The  Non-Computer  Design  Procedure 

1.  All  assumptions  for  the  prototype  shelter  e'ttend  to  the  design 
procedure. 

2.  The  air  temperature  transient  response  can  be  idealized  to  a  step 
function . 

3.  The  shape  of  the  shelter  does  not  seriously  affect  its  transient 
response . 
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Appendix  B 


CASE  2  METHOL  OF  ACCOUNTING  FOR 
MODEL-PROTOTYPE  TEMPERATURE  DISTORTIONS 


The  Case  2  method  Imposes  the  condition  that  here  is  a  one-to-one 
correspondence  in  heat  flux  between  the  model  'ind  prototype.  In  pursu¬ 
ing  this  approach,  both  the  model  shelter  wall  temperature  and  air 
temperature  will  be  distorted.  To  find  the  uist.>rtion  of  the  model  air 
temperature  when  h*  ft  consider  an  imaginary  model  for  which  h*  *  I, 
and  for  which  the  imaginary  wall  temperatur<i:  equals  the  actual  model 
wall  temperature.  The  imaginary  model  serves  as  an  intermediate  between 
the  actual  model  and  prototype,  and  it  yields  correct  air-wall  tempera¬ 
ture  differences. 


q 


P 


(B-1) 


h 

P 


(B-2) 


And 


t  =  t 
w  w 

m 


(t  -t  )  *  (t  -t  )  »  h*(t  -t  ) 

a  w  p  a  w  m.  a  w  m 


where:  q'  *a/A«h(t  -t  ) 

aw 


subscript  m^  refers  to  corrected  model 
temperatures  (Imaginary  model) 

subscript  m  refers  to  distorted  model 
temperatures  (actual  model) 

h*  =  h  /h  =  h  /h 

m  p  m  m^ 

From  Equation  B-4 


(B-3) 


(8-4) 


(t  -t  ) 
a  w' 


m. 


a 


m 


h* 


+  t 


w 


m 


t  -  t  +  h*t 
aw  w 

m.  m .  m 

i  t _ 

h* 


(B-5) 


(B-5a) 
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t  “  by  definition 
w  W 
m 

+  t  (h*-l) 
a  w 


i 

(B-5b) 

a  h* 

m 

h  (t  -  t  ) 

m.  a  w 

(B-6) 

1  m^  m 

.  t  - 

a_  h. 

(B-7) 

m  m. 


Substituting  Equation  B-7  into  B>5b  and  rearranging,  we  have 

t 


a 


m. 


t  =”T^“+ 
a  h*  ah 

m  m.  m, 

i  1 


-)  (1-  h?) 


(B-8) 


“i 

t  “  -TT—  +  t 
a  h*  a 
tn  n. 


m. 


h*  h  h  h* 
“i  “i 


^a  •  +  q'  (h  h*  ’  h  ^ 

IQ  IQ^ 


(B-8a) 


(B-8b) 


Now 


h  h* 
”i 


h  (r^)  -  h  -  h 
tn,  h  '  P  h  '  w 
ip  p 


(B-9) 


because  h  *  h  by  definition 
\  P 

Substituting  the  result  of  Equation  B-9  into  B-8b  and  rearranging,  we  have 

(B-8c) 


h  “h 
f-E — 5h 


t  -  t  +  q'  (r^^) 
a  a  ^  h  .h  ' 
mm  pm 


Note  that  the  last  term  of  Equation  B-8c  is  the  air  temperature  distortion, 

and  that  this  distortion  is  equal  to  zero  when  h  >■  h  . 

pm 


Remember  that  the  determination  of  the  prototype  air  temperature 
from  the  corrected  model  temperature  requires  that  the  distortion  of 
the  model  shelter  wall  temperature  be  found.  Here  again,  the  Imaginary 
model  concept  Is  useful. 


Now,  by  definition 


a 


+  (t 


'w>p 


<v 


t ) 

w  m. 


(t 


t  ) 
w  p 


(B-10) 


(B-11) 


Thus 


t 

a 


P 


t 

w 

P 


+  (t  -  t  )m 
a  w  1 


Since  q'*  *  1  for  the  Case  2  condition 

k  * 
w 

t  ®  i  t 
w  L*  w 
p  m 


(B-12) 


(B-13) 


(recall  that  t  “  t  ) 
w  w  ' 
m  m^^ 


Substituting  from  Equation  (B-13)  Into  (B-12)  and  rearranging 


k  * 

+  (jf  -  i>t 


w 

m 


(B-15) 


Now  the  Imaginary  model  air  temperature,  t  ,  can  be  eliminated  from 


Equation  (B-14)  by  substitution  from  (Equation  B-8c) 


a 


m 


h  -h  k  * 
pm  m 


(B-15) 


Since  the  model  may  not  have  the  same  Initial  temperature  as  the 
chosen  prototype,  all  temperatures  should  be  replaced  by  temperature 
differences  between  temperatures  at  time  0  and  the  Initial  temperature 
at  time  6  -  b. 


(t  -  t  ) 

Q  op 


h  -  h  k  * 

(t  -  t  )  -  q  *(..£,  ”)  +  (-—-  1)  (t  -  t  ) 

•e  ‘o  ''e  ''o 


(B- 
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Appendix  C 


BASIC  DATA  AND  CALCULATIONS  FOR 
SETTING  UP  THE  MODEL  RUN  --  WASHINGTON,  D.  C. 


The  basic  input  data  for  the  model  run,  Washington,  D.  C.  suitiner 
conditions  —  corresponds  to  the  values  used  by  Drucker  for  his  analog 
computer  solution  for  the  same  location.  This  procedure  allows  a  compari¬ 
son  of  the  model  behavior  to  the  analog  computer  solution. 

Temperatures  and  Absolute  Humidity 

Average  air  temperature  outside  shelter:  t^  ®  77°F 

Average  air  absolute  humidity  outside  w.  -  0.0162  Ib/lb 

shelter: 

Initial  soil  temperature(initial  air  t  *  SS^^F 

temperature  inside  shelter  equals  ^o 

initial  soil  temperature) : 


Air  Properties 
Specific  heat: 
Density: 

Soil  Properties 
Thermal  conductivity; 

Specific  heat: 

Density: 


c  »  0.241  Btu/lb-op 
a 

p  «  0.075  Ib/ft^ 
a 


k  =  0.83  Btu/hr-ft-op 
w 

P 

c  “  0.30  Btu/lb-Op 
w 

P 

p  -  140.0  Ib/ft^ 
w 


Other  Parameters 

Convective  heat  transfer  coefficient: 
Enthalpy  of  condensation: 

Total  shelter  area: 


h  -  1.47  Btu/hr-ft2-°F 
P 

H  -  1054  Btu/lb 

A  -  2884  ft^ 

P 
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Length  of  shelter  occupancy: 
Total  heat  output/person: 
Ventilation  rate/person 
Air  conditioning  rate/person 
Sensible  heat  factor: 

Number  of  persons  in  shelter: 


0  *10  days 

max 

*  400  Btu/hr 
g 

0*6  cfm 
R  -  0 
F  -  0.75 
n  >  100  people 


Model  parameters  were  equal  to  the  prototype  values  except  as  follows: 
Temperatures 

t 


Initial  soil  temperature  (Initial 
air  temperature  was  same  as  initial 
soil  temperature) : 

Soil  Properties 
Thermal  conductivity: 

Specific  heat: 

Density : 


Convective  heat  transfer  coefficient 
(determined  by  tests  with  constant 
heat  input) : 

Total  shelter  air: 


w 


72°F 


m 


w 


m 


w 


m 


m 


0.167  Btu/hr-ft-°F 
0.189  Btu/lb-OF 
92.3  Ib/ft^ 

1.15  Btu/hr- ft^-OF 


A  -  69.8  ft^ 
m 


At  this  point,  a  number  of  parameters  required  by  Equation  24b  can 
be  computed.  (All  computations  on  a  per  person  basis  where  applicable.) 

(hA)  ■  1.47  X  2884/100  "  42.4  Btu/hr-°F/person 
P 

L*  -  (A  /A  )^  -  (69.8/2884)^  -  0.155 
m  p 

h*  -  h  /h  -  1.15/1.47  -  0.782 
m  p 
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0.167/0.83  -  0.201 


k  * 
w 


k  /k 
w...  w 


h*L*  ,  0.782x0.155 

“  inf  “  0.201 


0.605 


( 


1-Bi*. 


.1-0.605. 
^  0.605  ^ 


0.654 


Assumliig  that  the  prototype  shelter  in  Washington,  D.  C.  will 
experience  an  air  temperature  rise  from  68°F  to  85°F,  the  heat  balances 
(Equations  25  and  26)  can  be  solved  for  various  dry-bulb  temperatures 
within  this  range.  The  values  obtained  thereby  are  the  prototype  heat 
inputs/person  as  a  function  of  prototype  air  temperature. 

’  ^®35-8t^)  -  R(F)  -  P,c^Q(W  (25) 

P 

q  -  (8t  -^35)  -  R(l-F)  -  P  0H(w  -  w  )  (26) 

a  p  a  a^^^  a  p 

when  q  *>0  (27) 

la 

P 

when  q^  «0  (28) 


Now  q^  is  the  total  heat  load  to  the  prototype  shelter  wall;  for  the 
P 

model,  the  heat  load  is 

qi  “  «V^*qX  (37) 

m  p 

-  100  X  0.201  X  0.155q^ 

^P 

■  3.13q^  Btu/hr 
P 


Recalling  that  (t  -t  )  ■  (t  -t  )  for  the  Case  1  condition,  Eviuation  15 

w  w  p  w  w  m  ‘ 

o  o 

requires  that: 


q*  ■  k  *L* 
^  w 


(C-1) 
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Since  q*  *  q  /q 


*  *^w*^*^p  ^  person  basis  (C-2) 


and  the  total  heat  input  to  the  model  is: 


qj  *  (C-3) 

m  p 


»  100  X  0.201  X  O.lSSq^ 


“  Btu/hr 

.  P 


*  0.919q^  Watts 
P 


Using  the  parameters  calculated  previously,  the  distorted  model 
temperature  rise  corresponding  to  each  prototype  temperature  can  be 
found  from  Equation  24b  rearranged  as 


9t  9>t<  1  _R'{4r 

<  V  ‘  ’p  * 

wo  Wo 


<  V  >n.  •  ( V  ‘O) "  <^>p  * 

0  0  0  ‘ 


where:  68  2 1  ^85 

“  “0  “ 


(24b) 


(C-4) 


The  several  solutions  of  equation  group  25,  26,  C-3,  and  C-4  for  q.,  q  , 

d  la 

q..  I  And  (t  -  t  )  can  be  conveniently  performed  on  an  IBM  1620  digital 
m  0  0 

computer.  Once  found,  the  distorted  model  temperatures  are  plotted  as  a 
function  of  q^  expressed  In  Watts  (see  Figure  3).  This  curve  Is  used 
m 

during  the  model  run  to  determine  the  sensible  heat  Inputs  associated 
with  distorted  model  air  temperatures.  Following  the  model  run,  the 
prototype  air  temperatures  are  found  from  the  distorted  model  air  tem¬ 
peratures  by  means  of  Equation  24b  (C-4).  It  Is  convenient  to  plot 
t  as  a  function  of  (t  -  t  )  as  shown  In  Figure  4. 

*0  *0  *0  ® 

P 
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Since  experimental  data  is  taken  in  real  time,  the  time  of  each 
temperature  observation  must  be  converted  to  prototype  time.  This 
conversion  is  made  from  Equation  10. 


e 

p 


a  * 
w 


a  * 


k  * 
w 


P  *c  * 


w 


0.201 

0. 923x0.66 
=  0.33 

=  A*  »  (69.8/2884) 

0.33x2884 

H  ®  — — — 

p  69.8 

=  13.6  6  (hours  for  each  hour  of  model  time) 
tn 

*  0.566  0  (days  for  each  hour  of  model  time) 


(10) 
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t,,  -  680p 
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« 


(°F) 
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figure  7.  Tanparature  Incramnts.  (I)  lacraoMats  are  added  to  corrected 
nodel  teaiparature  data  to  find  naxlauun  dry*bulb  temperature 
attained  during  the  day.  (2)  For  noQ->computer  design  procedure. 
Increments  arf  subtracted  from  t  allowed  to  find  average  air 

max 

temperature  attained  at  the  end  of  shelter  occupancy. 


FIG.  8  MODIFIED  PSYCHROMETRIC  CHART 


POUNDS  OF  WATER  PER  POUND  OF  DRY  AIR 
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(hours) 

riKuro  9.  Tronslont  r;osponso  of  tho  aodal  tholtor  wall  to  a  constant 
host  input. 


theoretical 


Figure  10.  Toaiparatura  dlatrlbutiana  la  the  eoll  surrounding  the  aodel  shelter  after 
8  hours  (conataot  heat  Input). 
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